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We use the instantaneous normal mode approach to pro- 
vide a description of the local curvature of the potential en- 
ergy surface of a model for water. We focus on the region of 
the phase diagram in which the dynamics may be described 
by the mode-coupling theory. We find, surprisingly, that the 
diffusion constant depends mainly on the fraction of direc- 
tions in configuration space connecting different local minima, 
supporting the conjecture that the dynamics are controlled by 
the geometric properties of configuration space. Furthermore, 
we find an unexpected relation between the number of basins 
accessed in equilibrium and the connectivity between them. 
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The study of the dynamics in supercooled liquids has 
received great interest in recent years due to possi- 
bilities opened by novel experimental techniques de- 
tailed theoretical predictions , and by the possibility of 
following the microscopic dynamics via computer simula- 
tion [Q. One important theoretical approach developed 
recently is the ideal mode coupling theory(MCT) 
which quantitatively predicts the time evolution of cor- 
relation functions and the dependence on temperature 
T of characteristic correlation times. Unfortunately, the 
temperature region in which MCT is able to make pre- 
dictions for the long time dynamics is limited to weakly 
supercooled states Q . In parallel with the development 
of MCT, theoretical work |^|-^ has called attention to 
thermodynamic approaches to the glass transition, and 
to the role of configurational entropy in the slowing down 
of dynamics. These theories, which build on ideas put 
forward some time ago |ic|-[l^, stress the relevance of 
the topology of the potential energy surface (PES) ex- 
plored in supercooled states. Detailed studies of the PES 
may provide insights on the slow dynamics of liquids, and 
new ideas for extending the present theories to the deep 
supercooling regime. 
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The instantaneous normal mode (INM) approach uses 
the eigenvectors (normal modes) and eigenvalues of the 
Hessian matrix (the second derivative of the potential en- 
ergy) to provide an "instantaneous" representation of the 
PES in the neighborhood of an equilibrium phase space 
point. Examples of INM analysis have appeared for a few 
representative liquids [^ - [l^ . In this Letter we present 
INM calculations for a rigid water molecule model, the 
SPC /E potential (l7| , for several temperatures and densi- 
ties in supercooled states. The reasons for selecting such 
a system are twofold: (i) The SPC/E model has been 
studied in detail, and it has been shown that its dynam- 
ics follows closely the predictions of MCT . Fur- 
thermore, the (p, T) dependence of the configurational 
entropy has been calculated, and shown to correlate with 
the dynamical behavior [Q. (ii) The SPC/E model has 
a maximum of the diffusion constant D under pressure 
along isotherms, as observed experimentally ||2^. The 
maxima in D can be used as sensitive probes for test- 
ing proposed relations between the change in topology of 
configuration space and the change in D. 

We analyze recent simulations for six densities be- 
tween 0.95 and 1.3 g/cm^ and five temperatures for each 
density. For each state point, we extract 100 equally- 
spaced configurations p^ , from which we calculate and 
diagonalize the Hessian using the center of mass and the 
three principal inertia momenta vectors as molecular co- 
ordinates j23|. For each configuration, we classify the 
imaginary modes into two categories: (i) shoulder modes 
and (ii) double-well modes, according to the shape of 
the PES along the corresponding eigenvector |^^. The 
shoulder modes are modes for which the negative cur- 
vature appears as a result of a local anharmonicity; the 
double-well modes include only the directions in config- 
uration space which connect two adjacent minima. The 
classification of the modes is done by studying the shape 
of the PES along the eigenvector corresponding to the 
mode, i.e. beyond the point in configuration space where 
such modes have been calculated. In this respect, the 
potential energy along the eigenvector direction may be 
different from the actual curvilinear direction described 
by the time evolution of the eigenvector under considera- 
tion. Furthermore, some modes classified as double-well 
may be related to intra-basin motions | |25[ | and therefore 
may not be relevant to describe diffusivity. In the present 
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study we show — using the D maxima hne as a guide — 
that these spurious features are neghgible. 

Fig. |l] shows the density dependence of the fractions of 
imaginary (unstable) modes /„, double- well modes fdw, 
and of shoulder modes fsh, where = fdw + fsh] we also 
show D for the same state points. The close correlation of 
fdw and D is striking. At p k, 1.15 g/cm'^, D has a max- 
imum; at this p, the increase of diffusion constant caused 
by the the progressive disruption of the hydrogen bond 
network on increasing the density is balanced by slowing 
down of the dynamics due to increased packing. At the 
same density, /„ and fdw also show maxima, supporting 
the view that these quantities are a good indicator of 
the molecular mobility. There is also a weak maximum 
in fsh, but at a p < 1.15 g/cm^ j2^. The presence of 
maxima in both D and fdw at the same density suggests 
that, for the SPC/E potential, fdw is directly related to 
D. The relation between D and fdw is shown in Fig. || 
for all the studied isochores. We note that D in this sys- 
tem is a monotonic function of fdw and that all points 
fall on the same master curve [|7|, notwithstanding the 
large range of D values analysed. Thus, surprisingly, in 
SPC/E water, the knowledge of the fraction of directions 
in configuration space leading to a different basin, fdw, is 
sufficient to determine the value of D. We also note 
from Fig. ^ that D vanishes, but at a small non-zero value 
of fdw ~ 0.007 (vertical arrow) ||2^, suggesting that in- 
deed a small number of spurious double-well modes, re- 
lated to intrabasin motion [^,^, are still included in 
our classification. 

The reduction of mobility on cooling in the studied 
(p, T) range - where MCT provides a description of the 
dynamics - appears related to the geometrical properties 
of the PES, i.e. the system mobility is reduced because 
the number of directions connecting different local min- 
ima (needed to explore freely the configuration space) is 
decreasing [ ^lljs^ ]. Hence, the observed reduced mobility 
is mainly related to the geometry of the PES, i.e. it is 
"entropic" in origin. 

MCT seems to be able to describe the entropic slow- 
ing down of the dynamics associated with the vanishing 
of fdw- This is consistent with the general consensus that 
the missing decay channels for the correlations responsi- 
ble for the failure of MCT at very low T are activated 
processes. To make closer contact with MCT, we com- 
pare in Fig. ^the density dependence of the MCT critical 
temperature Tmct jisl with the T at which the fraction 
of double-well modes goes to zero, as well as with the T 
at which fdw reaches the estimated asymptotic value of 
0.007 ||3^, to take into account the small over counting 
intrinsic in the mode classification. We observe that the 
fdw = locus tracks closely the Tmct hne, but that the 
fdw — 0.007 locus nearly coincides with the Tmct hne. 

The results presented here suggest that the liquid dy- 
namics in the MCT region of SPC/E water is controlled 
by the average connectivity in configuration space. On 



the other hand, for the same model, the (p, T) depen- 
dence of the configurational entropy .Sconf ~ which, in the 
Stillinger- Weber |l^ formalism, can be defined as the 
logarithm of the number of different basins f2 in config- 
uration space - has also been calculated and shown to 
correlate with D ]20| . Since no a priori relation is ex- 
pected between the connectivity of the local minima and 
their number, we consider particularly interesting the ob- 
servation of such a relation. For this reason, we show in 
Fig. ^ a "parametric plot" (in the parameter D) of fdw 
versus ^conf for all studied densities. Fig. ^ shows a lin- 
ear relation between ln{fdw) and Sconi, i-e. between fdw 
and the number of basins (since = exp(S'conf/fcB)) and 
highlight the existence of an unexpected relation between 
the number of basins accessed in equilibrium and the con- 
nectivity between them. This novel feature of the PES 
contributes to a better understanding of the water dy- 
namics, and may contribute to the understanding of the 
dynamics of glass- forming liquids [^3| . 
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FIG. 1. Density dependence along the studied isotlierms of 
(a) fu, (b) fdm, (c) fsh- Temperatures are T = 21QK (circles), 
T = 22QK (squares), T = 22.QK (diamonds) and T = 240^" 
(triangles) ._(d) shows the p dependence of D for the same 
isotherms ]19[. 



10" 



^10 

E 
Q 

10 



10" 




p = 0.95 g/cm 
p = 1.00 g/cm 
p = 1.05 g/cm 
p = 1.10 g/cm 
p = 1.20 g/cm^ 
p = 1.30 g/cm 



0.00 



0.01 



0.02 



0.03 



dw 



FIG. 2. Diffusion constant D versus fdw for different iso- 
chores. The arrow indicates the value fdw ~ 0.007, to high- 
light the presence — even when D approaches zero — of a 
small number of spurious double- well modes, related to intra- 
basin motion paBOl. 
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FIG. 3. Comparison of the mode coupling critical temper- 
ature Tmct [M , of the fdw ~ locus and of the fdw = 0.007 
locus. 
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